Non Fermi Liquid Dynamics of the Two-Channel Kondo Lattice 
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The paramagnetic phase of the two-channel Kondo lattice model is examined with a Quantum 
Monte Carlo simulation in the limit of infinite dimensions. We find non- Fermi-liquid behavior at low 
temperatures including a finite low-temperature single-particle scattering rate, no Fermi distribution 
discontinuity, and zero Drude weight. Both the optical and quasiparticle mass enhancement and 
scattering relaxation rate show consistent evidence of non-fermi liquid behavior. However, the low- 
energy density of electronic states is finite. 
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Introduction. The Fermi Liquid theory of Landau has 
provided a remarkably robust paradigm for describing 
the properties of interacting Fcrmion systems such as liq- 
uid 3 He and alkali metals (e.g., sodium). The key notion 
of this theory is that the low lying excitations of the inter- 
acting system possess a 1:1 map to those of the noninter- 
acting system and hence are called "quasiparticles." In 
the metallic context, one may think of the quasiparticles 
as propagating electron-like wave packets with renormal- 
ized magnetic moment and effective mass reflecting the 
"molecular field" of the surrounding medium. A sharp 
Fermi surface remains in the electron occupancy func- 
tion nj: which measures the number of electrons with a 
given momentum, and for energies u> and temperatures 
T asymptotically close to the Fermi surface the excita- 



tions have a decay rate going as u> 



! (k B T) 2 , which 



is much smaller than the quasiparticle energy. This gen- 
erally translates into a T 2 contribution to the electrical 
resistivity p(T) and an lu 2 contribution to the scattering 
relaxation rate T(uj) inferred from the optical conduc- 
tivity, and both the optical and single-particle electronic 
masses are found to be enhanced but positive and finite. 

The Fermi liquid paradigm appears now to be break- 
ing down empirically in numerous materials, notably a 
number of fully three dimensional heavy Fermion alloys 
and compounds In these systems such anomalies as 
a conductivity with linear dependence on uj, T and loga- 
rithmically divergent linear specific heat coefficients are 
often observed. In addition, some heavy fermion lattice 
compounds, such UBei3 Q and CeCu2Si2 j|, the super- 
conductivity arises in a normal state which is clearly not 
described as a Fermi liquid. Specifically, the resistivity is 
approximately linear in T above the transition, and the 
magnitude of the resistivity at the transition is high (typ- 
ically 80-100 jti— 51-cm in UBei3 even in the best samples). 
If the quasiparticle paradigm indeed breaks down, this 
may require completely new concepts to explain the su- 
perconducting phases of these materials. While the Lut- 
tinger liquid theory provides an elegant way to achieve 
non- Fermi liquid theory in one-dimension (with, e.g., no 
jump discontinuity in nr, and separation or unbinding 



of spin and charge quantum numbers), this results from 
the special point character of the Fermi surface. Whether 
the essential spin-charge separation may "bootstrap" into 
higher dimensions remains unclear Among the re- 
maining theories to explain experiment are those based 
upon proximity to a zero temperature quantum critical 
point |^| , those based upon disorder induced distributions 
of Kondo scales in local moment systems , and those 
which hope to explain the physics from impurity to lat- 
tice crossover effects in the multi-channel Kondo model 
[Q. Notably lacking for dimensions higher than one are 
rigorous solutions to microscopic models which display 
non- Fermi liquid behavior . 

In this paper, we present the first rigorous solution for 
the dynamics of the two-channel Kondo lattice model in 
infinite spatial dimensions. We find that the paramag- 
netic phase of this model is an "incoherent metal" with 
finite density of states at the Fermi energy and finite 
residual resistivity. Both the single-particle and the op- 
tical excitations are non-Fermi liquid like; in particular, 
there is a finite lifetime for electrons at the Fermi energy, 
an ill defined quasiparticle mass, a linear low tempera- 
ture electrical resistivity with a finite residual value, a 
negative optical mass and a linear in frequency low tem- 
perature scattering relaxation rate. 

Motivation The two-channel Kondo impurity model 
consists of two identical species of non-interacting elec- 
trons antiferromagnetically coupled to a spin 1/2 impu- 
rity. Non-Fermi liquid behavior results because of the 
inability to screen out the impurity spin: it is energet- 
ically favorable for both conduction electron bands to 
couple to the impurity which gives a spin 1/2 complex 
on all length scales as the temperature goes to zero. As a 
result, the ground state is degenerate and the excitation 
spectrum non-Fermi liquid like. We believe that at least 
the quantitative features of the model will change once 
the concentration of impurities becomes finite. A single 
impurity always lives in a Fermi liquid host; whereas the 
introduction of a finite concentration of impurities into 
the system will mean that at some temperature, the host 
of each impurity will be a non Fermi liquid due to the 
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presence of the others. The simplest model to test this 
hypothesis is the fully concentrated one; ie. the two- 
channel Kondo lattice model. 

Model The Hamiltonian for the two-channel Kondo 
lattice is 



z.a.cr 



(1) 



where cj Q CT [ci ia>a ) creates (destroys) an electron on site 
i in channel a = 1, 2 of spin cr, S,; is the Kondo spin on 
site i, and Sj )Ct are the conduction electron spin operators 
for site i and channel a. The sites i form an infinite- 
dimensional hypercubic lattice. Hopping is limited to 
nearest neighbors with hopping integral t = t* /2y/d; the 
scaled hopping integral t* determines the energy unit and 
is set equal to one (t* — 1). Thus, on each site the Kondo 
spin mediates spin interaction between the two different 
channels. This problem is non trivial, and for the re- 
gion of interest in which J > and T « J, C it is de- 
scribable only with non-perturbative approaches. Clearly 
some simplifying method which allows for a solution of 
the lattice problem in a non-trivial limit is necessary. 

Formalism and Simulation Such a method was pro- 
posed by Metzner and Vollhardt || who observed that 
the renormalizations due to local two-particle interac- 
tions become purely local as the coordination number 
of the lattice increases. A consequence is that the so- 
lution of most standard lattice models may be mapped 
onto the solution of a local correlated system coupled to 
an effective bath that is self-consistently determined [|l0| . 
We refer the reader to the above references and recent re- 
views for further details on the method jy]]. In order to 
solve the remaining impurity problem, we use the Kondo 
impurity algorithm of Fye and Hirsch flSql , modified to 
simulate the two-channel problem Jl3| . The results pre- 
sented here are limited to the model at half filling of the 
conduction bands (JV = 1.0 for J = 0.75,0.625,0.5,0.4). 
A sign problem was also encountered in the QMC process 
which limited our access to very low temperatures. The 
Euclidean-time QMC results for the local greens function 
G(t) were then analytically continued to real frequencies 
using the "annealing" Maximum Entropy method [p"4[ . 
The single particle self energy may then be obtained by 
inverting the relation G(u) = —i^/irw (uj + fi — , 
where w(z) is the complex Fadeev function. Error bars 
are ill-defined for individual points in analytically con- 
tinued spectra O], and are less than 3% for the other 
results presented here. 

Results. When the temperature is lowered below the 
lattice Kondo temperature T Ea], we find non- Fermi 
liquid behavior in the single-particle properties of the 
model. For example, the derivative of the particle dis- 
tribution function dn(ek)/dek = —T^ n l/ (iuj n — + 



fi — S(iw„)) 2 saturates to a finite width distribution p5| , 
in contrast to a Fermi liquid where it would display a 
dominant delta function contribution as T — > 0, and a 
Luttinger liquid Q or Marginal Fermi liquid fl6[ | where 
it would have a singular divergence. 

In Fig. 1 we show one electron properties of the model. 
Fig. 1(a) displays the single particle DOS. This has a 
finite value as lu,T — > 0, with a peak away from the 
Fermi energy. Novel behavior is seen in the real part 
of the one electron self energy which has positive slope 
at uj — > (Fig. 1(c)). For a Fermi liquid, this slope 
would be negative. The physical content is important: 
Z = 1/(1 — dReTi/duj) measures the overlap of the quasi- 
particle wave function with the original one-electron wave 
function having the same quantum labels. A positive 
slope leads to Z > 1 or Z < 0, indicating a breakdown 
of the quasiparticle concept. Concomitant with the fi- 
nite width of —dn/dek is a finite imaginary part to the 
low temperature self energy (Fig. 1(b)). This indicates 
that the one-electron excitations arc ill defined on ap- 
proach to the Fermi surface, again ruling out a Fermi 
liquid description. Since the low temperature thermody- 
namic properties such as the specific heat, uniform mag- 
netic susceptibility, and charge susceptibility display no 
evidence for a gap, we believe the observed behavior in- 
dicate a new kind of metallic state. 

From our numerical results, one can show that the zero 
temperature self energy must be non-analytic. If the self 
energy is analytic everywhere in the upper half complex 
plane (and infinitesimal close to the real axis), then one 
may easily show that 



lim = lim ; 

T^0 LU T-+0 duj 



(2) 



where lu = ttT is the lowest Fermionic Matsubara fre- 
quency along the imaginary axis. However, we find that 
when T < To, Im ^"° - < and is apparently divergent; 

whereas, from Fig. 2 it is apparent that 



dR.cSf 



duj 

1.5. Thus, S(w) cannot be analytic at the origin of the 
complex plane. 

The non-Fermi liquid behavior also strongly effects ex- 
perimentally relevant transport properties of the system. 
The electrical resistivity is shown in Fig. 2. We find that 
p(T)/p(0) curves for different J values collapse onto a 
universal scaling curve when plotted against T/Tq. As 
shown in the inset, p(T) ~ p(0)[l + B(T/T )} for T -> 0, 
with B < 0. We interpret the finite value of p(0) to- 
gether with — im£(0, 0) as "spin disorder scattering" off 
of the degenerate screening clouds centered about each 
local moment spin. Given the finite density of one par- 
ticle excitations at the Fermi energy, this finite residual 
resistivity is indicative of an "incoherent metal phase" 
brought about by the disordered spin degrees of free- 
dom, in qualitative agreement with results obtained with 
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a Lorentzian bare conduction DOS (which doesn't self- 
consistently renormalize) Jl7| . This surprising fact of a 
finite resistivity for T — > in a translational invariant 
system reflects the infinitely degenerate ground state of 
the two channel lattice with a residual entropy. 

The optical conductivity ui(lu) is measured in 
units of (To = e 2 ir/2ha, which varies between 
10~ 3 ...10 _2 [(/ificm) _1 ], depending on the lattice con- 
stant a. As shown in Fig. 3(a), consistent with the 
lack of quasiparticles, the low temperature optical con- 
ductivity displays vanishing Drude weight D. To con- 
firm this, we also calculated D by extrapolation [ fL8| , 
D oc Xjiio. = °> iuj a ->■ 0) - Xjj{l = 0, iuJ a = 0), where 
Xjj(q = 0,itJ a ) is the bulk current-current susceptibility. 
For each value of J and all fillings, we found that the ex- 
trapolated zero-temperature Drude weight was zero. In 
addition, ai(u>) has a finite-frequency (u> 0.6 J) peak. 
Both these features again support our interpretation in 
terms of a new kind of non-Fermi liquid metallic state. 

The optical conductivity of metals, even non- Fermi liq- 
uid metals |pp| , is usually analyzed by rewriting it in a 
generalized Drude form 



cr(w) 



f 



(3) 



where <j(lj) = ai(u>) + ia%{u)). The resulting mass 
enhancement (1 + A(w)) and scattering relaxation rate 
r(ct>) are shown in Fig. 4 when J = 0.625 for various 
temperatures. The lifetime approximation becomes ex- 
act in infinite dimensions, thus at low u> one expects 
1 + X(uj) = 1 - ^Re£(w/2) and T{uj) = -2ImS(w/2) 
p9| . At high temperatures, and (1 + A(w)) are 

essentially flat so a Drude peak is recovered in ai(u>); 
whereas at lower temperatures they become strongly fre- 
quency dependent. Here, the behavior of T(uj) is con- 
sistent with the imaginary part of the self energy, indi- 
cating that the optical and quasiparticle scattering rates 
are consistent for this model, even though quasiparticles 
do not exist. Note that when T/To = 0.20, the low fre- 
quency optical mass is negative, m*(0) » —1.0. This 
behavior is only qualitatively consistent with the mass 
enhancement estimated from the low temperature self 
energy 1 - dReS(w)/dw| w=0 w -0.5 (cf. Fig. 1). For 
larger values of J > 0.75, m*(0) > 0. Thus, the negative 
mass enhancement should be viewed as a sufficient but 
not necessary identifying feature of a non Fermi Liquid. 

Finally, we mention the possible applicability of our 
results to concentrated heavy Fermion systems. Three 
systems display resistivity of the form p(T) w p(0)[l + 
B(T/T ) a ] for T < T , with p(0) of order the unitarity 
limit and a ~ 1. They are: UCus-^Pd^ (with B < 0) 
[|, UBei 3 (B > 0) |H,||, and CeCu 2 Si 2 (B > 0) 
[3|. UBei3 and CeCu2Si2 are ordered compounds which 
have been proposed as possible two-channel lattice sys- 
tems (see Refs. @(b,c)), while UCu5_ x Pd x is a possible 



example of the distribution of Kondo scales scenario || . 
In both UBei3 Q and UCu 5 _ a ,Pd x |Q the measured 
<7i(lu) and T(u>) are strikingly similar to those displayed 
in Figs. 3 and 4. 
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FIG. 2. Resistivity of the two-channel Kondo lattice. 
p(T)/p(0) is plotted versus T/To for various values of J. In 
the inset, the lowest temperature data (for T/Tb < 1, shown 
as open circles) was fit to p(T)/p(0) = 1 + B(T/T ) a , with 
B = -0.4 and a = 1.03. 



S 




CD 




2.0 


T/T =0.20 


O 


1.5 


T/T =0.41 

T/T =0.82 

T/T =1.63 


3 

0~ 


1.0 
0.5 


T/T =3.27 




0.0 




FIG. 1. Single particle properties of the two-channel Kondo 
lattice in infinite dimensions when J = 0.625 and N = 1.0. (a) 
Single-particle density of states (DOS). At high temperatures 
T To (not shown), the DOS is a Gaussian, crossing over 
to the peaked distribution with relative suppression at u — » 
for lower temperatures when T -C T). (b) Imaginary part of 
the self energy. As the temperature is lowered the self energy 
does not approach a Fermi liquid form ImE(u) oc — T 2 — lo 2 , 
but rather appears to approach the non-analytic form (see 
text) ImE(oj) oc — c + (c) The real part of the self en- 
ergy ReE(oj), is also anomalous since its initial slope is posi- 
tive indicating a quasiparticle renormalization factor which is 
greater than one. 
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FIG. 3. Optical conductivity of the two-channel Kondo lat- 
tice, (a) Optical conductivity when J = 0.625. As the tem- 
perature is lowered T < To, the optical conductivity devel- 
ops a pseudogap at low frequencies. No evidence of a finite 
Drude weight, D, can be seen, consistent with estimates of D 
obtained from extrapolation of Matsubara-frequency results 
@. (b) The optical conductivity vs. uj/J when T = 0.0156 
for various values of J. a\(uj) displays a finite-frequency peak 
at roughly lo w 0.6J. 
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FIG. 4. Optical scattering rate r(w)(a) and mass enhance- 
ment (1 + (b) obtained from the optical conductivity 
shown in Fig. 2(a) using Eq. |. For T/T « 1 both T(u) and 
(1 + A(w)) display non fermi liquid behavior: r(w) is roughly 
linear in ui and m*(0) < 0. Typical Fermi liquid behavior 
is recovered when T/To ^ 1, since both F(ui) and (1 + A(w)) 
become weakly frequency dependent consistent with the re- 
covery of the Drude peak in Fig. 3(a). 
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